Light curve and neutrino spectrum emitted during the 
collapse of a nonrotating, supermassive star. 
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p • Abstract 

' The formation of a neutrino pulse emitted during the relativistic collapse of a 

spherical supermassive star is considered. The free collapse of a body with uniform 
density in the absence of rotation and with the free escape of the emitted neutrinos 
can be solved analytically by quadrature. The light curve of the collapsing star and 
the spectrum of the emitted neutrinos at various times are calculated. 

1 INTRODUCTION 

The properties of a neutrino pulse arising during the spherically symmetrical collapse of a 
^ , supermassive star into a black hole have been considered in a number of studies, the earliest 
' of which were carried out nearly 40 years ago [1]. It was shown in [2] that, in general, the 
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asymptotic time dependence of the total luminosity of the star has the form 
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; The propagation of photons outside a nonrotating collapsing star was considered in detail 
. ^ ' in [3] . The brightness distribution across the visible disk and the time dependence of the 
^ stellar spectrum were obtained for late stages of the collapse. The results of numerical com- 
putations of the light curve and the spectrum of a nonrotating, collapsing star are presented 
by Shapiro [HIH], who computed the propagation of photons inside the star in a diffusion ap- 
proximation taking into account general relativistic effects. The propagation of photons to a 
distant observer is also treated in detail in [3]. Estimates of the parameters of the spectrum 
and of the total intensity of a neutrino pulse arising during the collapse of a supermassive 
state were obtained in [6|. The possibility of detecting such a pulse using modern detectors 
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was analyzed, and the general background of neutrinos arising as a result of explosions of 
supermassive stars was estimated. The results of relativistic hydrodynamical computations 
of the collapse of a nonrotating, supermassive star are presented in [7], together with tem- 
perature and density profiles for the matter inside the star at various times. However, the 
propagation of the neutrinos was computed using a very simple scheme; it was essentially 
assumed that the neutrinos propagated instantaneously. In addition, all general relativistic 
effects were neglected. 

Thus, in spite of the abundance of publications on the topic, the properties of neutrino 
pulses arising during the collapse of supermassive stars are not fully understood. As a rule, 
previous papers have either considered very simple models for this phenomenon that have 
admitted exact solutions (as in [1]), constructed very detailed models that have required 
the use of approximation computational methods, or have neglected certain physically im- 
portant effects. For example, the propagation of the neutrinos is taken to be instantaneous 
in [7], and gravitational forces are not taken into account. Such an approach is clearly not 
justified physically. During the collapse of a supermassive star, the temperature of the mat- 
ter grows as the collapse proceeds; accordingly, the neutrinoloss rate in a comoving frame 
also grows, reaching its highest values in the late stages of the collapse. The drop in the 
radiation intensity at the end of the collapse is due exclusively to gravitational effects: the 
gravitational redshift and time dilation. Thus, the shape of the light curve is determined by 
the battle between two competing processes: the increase in the brightness of the star due 
to heating and the decrease in its brightness due to gravitational effects. In addition, the 

T 

characteristic time scale for the collapse — is comparable to the time for a neutrino to travel 

c 

a distance equal to the radius of the star. Under such conditions, the neglect of gravitational 
effects and of the neutrino-propagation time is completely unjustified. Ordinary (not su- 
permassive) stars were considered in [5]; these computations included all general relativistic 
effects, but the propagation of the neutrinos inside the star was treated in an Eddington 
approximation, which cannot be formulated unambiguously inside the star, where the effects 
of sphericity are appreciable. 

The aim of the current paper is to construct a model for the neutrino emission radiated 
by a supermassive star that can yield as accurate a solution as possible, including all general 
relativistic effects, while simultaneously providing a good description of a real collapsing star. 
The light curve and spectra obtained with this model should display all the main properties 
characteristic of the light curves and spectra of real collapsing supermassive stars, while the 
simplicity of the model enables us to understand the dependence of the parameters of the 
neutrino impulse on the parameters of the problem. As we will show below, by introducing 
a number of natural constraints on the parameters of the problem, we can obtain a solution 
by quadrature. 

2 DYNAMICS OF THE COLLAPSE 

We introduce a number of assumptions about the physical parameters of the collapse that 
simplify the computations. Most importantly, the star is taken to be nonrotating and the 
distribution of matter to be uniform at all times. In addition, we neglect completely the 
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influence of tlie pressure on tlie motion of tlie matter; i.e., tlie dynamics of the collapse are 
taken to be the same as for a sphere of dust. We also assume that the matter is heated 
according to a power law as a result of compression. Our neglect of the pressure is a rather 
crude approximation at the onset of the collapse, but the temperature of the matter is still 
modest at this stage, and few neutrinos are radiated. As the star approaches its gravitational 
radius r^, the influence of the pressure tends to zero, so that it indeed becomes negligible. 
As is shown in [8], inside a homogeneous dust sphere that is shrinking or expanding, it is 
possible to introduce a comoving, synchronous coordinate system in which ordinary three- 
dimensional space has the same curvature everywhere, which depends only on time. The 
character of this curvature (positive, negative, zero) does not change in the course of the 
collapse and is uniquely specified by the initial conditions (initial density and velocity). 
We will take the three-dimensional space inside the sphere to have zero curvature, which 
appreciably simplifies the computations without introducing large errors. 

Let the initial radius of the star be x^. According to [S], the metric inside the dust 
sphere can be written as 



ds^ = dr^ - a''{T)[dR^ + R^{dC + sin^C df)] (1) 
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where a(r) = f 1 - — ) ' (2) 

V Tcol J 

Here Tcou is a parameter of the problem that is defined below and the "spatial" coordinates 
{RX,0 ^^^6 analogous to the usual spherical coordinates in threedimensional space. T he 
coefficient a plays the role of a scale factor: to find the physical distance between two points 
in the metric ([1]), we must multiply the "coordinate distance" [calculated in the coordinates 
(i?, by a. In particular, the distance from a point with radial coordinate R to the 
center of the star is: 

X = aR (3) 

The coordinate system {t,R,(,0 is comoving. This means that each element of the col- 
lapsing matter has some fixed R, ( and ^ that do not change in time. In this sense, {R, CO 
can be considered Lagrange coordinates. 

Note the full analogy between the solution for a collapsing dust sphere considered here 
and the Friedmann solution for a fiat universe filled with homogeneous dust, whose metric 
coincides fully with[Tl Therefore, the inner space of the dust sphere can be considered part 
of a Friedmann universe; the motion of the neutrinos within the dust sphere occurs precisely 
as within a Friedmann universe. In a comoving coordinate system, the radial coordinate Rg 
of the surface does not change and at any time r < Tcoi is equal to its initial value 

Rs = — = Rso = -— (4) 
a ao 

In addition, it is straightforward to obtain the following relation using the results of [8]: 

27ike = (5) 



^Here and below, unless explicitly stated otherwise, we will work with a physical system of units in which 
c—1 and Vg = — 1, which specifies the units of velocity (c), length (r^), and time (rg/c). 
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Here, k is the gravitational constant and e is the density of the stellar matter in its rest 
frame. Integrating e over the entire volume of the collapsing sphere yields the total rest mass 
of the stellar matter. It is of interest to determine how this quantity is related to the total 
mass of the star measured by a distant stationary observer. This latter mass consists of 
the rest mass of the matter and the sum of its kinetic and potential energy. As was noted 
above, we are considering the collapse of a dust sphere within which the three-dimensional 
curvature of space is zero. In this case, the motion of dust matter toward the center is a 
free fall with zero velocity at infinity. The sum of the kinetic and potential energy of matter 
undergoing such a free fall is zero. Therefore, the mass of the star is equal to the rest mass 
of its matter, i.e., 

M = JedV (6) 
Hence, taking into account also the homogeneity of the star, we have 

3M 



Combining this relation with ([5]) yields 



(7) 



3kM 1 2R^ 

— = — 2 — other words, in accordance with (jlj) kM = — ^. (8) 

Sr^Q^a ^^coi 

Since in our units = c = 1, it follows from the definition Vg = that 

kM = l (9) 



which yields 



4/?3 2 3 

^ = 1, r^oi = -m (10) 



We will adopt the time r = as the initial time. Then, oq = 1 and 

XsO = RsO = Rs- (11) 

Having specified the initial radius of the star Xsq, we can use ( ITU]) to determine the collapse 
parameter Tcoi- The collapse of a uniform sphere with a flat internal space involves two 
independent parameters: the Lagrange radius of the surface Rg and the mass of the star. 
The quantity Tcoi — tq is uniquely determined by these two parameters, and variations in 
To only shift the zero time. In the case of arbitrary three-dimensional curvature, we can 
arbitrarily specify the mass, initial radius, and initial collapse velocity of the star. If we 
consider a collapse with zero curvature, this is not possible, and the initial collapse velocity 
always differs from zero for a finite initial radius. This is automatically determined from 
the condition of zero curvature and turns out to be equal to the velocity for free fall from 
infinity for the radius Xso- In the interest of conciseness, we will denote the initial radius of 
the star simply as Xq. 
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3 COMPUTATION OF THE RADIATION SPECTRUM 
IN A COMOVING FRAME 

We will assume that each matter element emits neutrinos isotropically in a comoving frame 
and that the intensity and spectrum of the emitted neutrinos depend on the temperature 
and density of the matter. In our analysis, the dynamics of the collapse corresponds to the 
motion of the dustlike material. If the law for variations of the density and the parameters 
M and Rs is known, the variations of the temperature in the compressed matter can be 
found from the solution of the energy equation for a specified initial entropy using the known 
rate of neutrino emission as a function of the temperature and density. Below, we use a 
simplified approach in order to obtain an analytical solution, in which we approximate the 
time dependence of the temperature with a power-law dependence of the temperature on 
the density. The variations of the density of the matter e(r) are determined in accordance 
with (jl]) and ([7]) by the single function a(r). To find the observed light curve and spectrum 
of the neutrinos, we must first determine the propagation toward the observer of a neutrino 
signal arising within a uniform collapsing body of a specified mass M. 
We will suppose that 

f{a,q)dV£pdT (12) 

neutrinos are emitted in an element of physical^ volume dV in an element of physical phase 
space d^p in an interval of proper time dr in the element's rest frame, where /(a, q) is some 
function whose form will be determined below. 

We introduce the neutrino distributioijfl function N such that the quantity Nd^pdV 
specifies the number of neutrinos in an element of physical volume dV in an element of 
physical phase space d^p. In the spherically symmetric case with isotropic neutrino emission 
in the comoving frame, N depends on only four quantities: the radial coordinate R, the 
angle between the neutrino trajectory and the direction from the center of the star the 
neutrino energy q, and time r. In this cas^, the kinetic equation can be written in the form 

m 

dN sin^dN 2q dN 2JEdN ^, ^ 

^^^^IrB B-17T + ^ r~a q 1^ = «/(«,9) (13) 

oR R ov oTcowJa oq 6Tcoi oa 

Here, the parameter a, which depends on time in accordance with ([2]), is used in place of 
the time r. The characteristic system for this equation has the form 



dR 

cos-d 
_Rd'd 

'iTcoi\/adq 
2q 



dy (14) 
dy (15) 
dy (16) 



^Not coordinate. 

•^We will call a quantity defined in this way a physical distribution function. 

^We also apply here the condition of zero three-dimensional curvature of the space inside the star. 
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^^-'^^ dy (17) 



2v^ 

The solutions of this system determine the trajectories of individual neutrinos. Combining 
pairwise the equations for characteristic system flTB]) and ffT7|) . and also ffUl) and ffT^ . we 
obtain the two first integrals: 

A = aq (18) 
B = Rsm'd (19) 

Combining (ITBI) and (fT7|l and using (IT^ . we obtain the third, last, integral of the system 

(ED-dZD- 

C = 3noiVa + R cos ^ (20) 

The physical meaning of these quantities can be elucidated based on the analogy between 
the solution for a collapsing dust sphere and the solution for a Friedmann universe. As 
was noted above, the motion of a neutrino is absolutely identical in these two cases; in 
particular, all three of the first integrals of system f[T^ - fll7l) are also valid in a Friedmann 
universe. The presence of integral f|T8|) . i.e., the fact that the product a ■ g is constant, leads 
to the appearance of a "redshift" in a Friedmann universe. As the parameter a increases 
(corresponding to the expansion of space), the energy of a neutrino along its trajectory 
decreases. In contrast, in our case we have a compression of a sphere, so that a decreases 
and the neutrino energy increases. Thus, we can think of our case as giving rise to a "violet 
shift". 

Supermassive stars with masses exceeding ~ lO^M© have small optical depths to the 
absorption of neutrinos, which can be taken to propagate freely. Taking the radius of the 
star to be i? = nRg = 2nGM/c^, we obtain expressions for the mean density, 

3M 1.8 ■ 10^6 /M0\ , , 

g/cm 



and the optical depth to neutrino absorption, 



r. = 10-^Vi?=^^«l 
M 



for n ~ 3, M = 1O^M0. 

Let us elucidate the qualitative form of the characteristics determining system (|nil - (|T711 . 
The motion of photons and neutrinos is linear in a Friedmann universe. Indeed, integral 
( 1191) shows that a neutrino moves along straight lines in the (r, i?, C, coordinate frame. In 
addition, since this frame is comoving, the collapsing star can be represented at any time 
in this frame as a sphere with a constant radius Rg. The trajectory described by ( fT4l) -( fT7l) 
proves to be a chord of this sphere. In accordance with the definition of the angle 9, the 
trajectory enters the sphere at some angle '& > \ and leaves at an angle '(9 < f . Proceeding 
from this property of the characteristic, it is straightforward to impose the appropriate 
boundary conditions for (fT3!) . For this, it is sufficient to specify the value of the function 
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at the surface of the star for angles "i? > |. We will assume that no neutrinos impact the 
star from outside, in which case 

mR=R. = > I (21) 

This is the boundary condition for f[T^ . 

We will denote the initial point of the characteristic and all related quantities with the 
subscript I and the final point and related quantities with the subscript II. Since both points 
are located on the surface of the collapsing sphere, 

Ri = Rii = Rs (22) 

Hence, we obtain using (fT9|) : 

sin 'di = sin 'djj (23) 

Since > f,^?// < f, 

cos^i = — cos^ii (24) 
In addition, it follows from the existence of the integral ( ITSl) that 



aiqi = aiiqii .. aj = an— (25) 

Finally, writing integral (120]) . we obtain 

3Tcozya7 + Ri cost?/ = STcoiy^on + Rii cos^u (26) 
Using relations f l2^ and ( l25i) . we obtain 



\aii— = Vau + ^^cos^n (27) 
V qi ^Tcoi 

Qi = — ^ (28) 

2R, 



cos-i? 



II 



\ ^'^coly^oJi 

Let us turn to the solution of (|T3l) itself. We first take the neutrino-emission function f|T2l) 
to be 

f{a,q) = Xa^5{q~qo) (29) 

i.e., each element of matter in the proper frame radiates monochromatic neutrinos with 
energy go- Here, A and (3 are numerical parameters. Since ( |T3l) is linear, the resulting 
solution can easily be generalized to the case of a more complex function f{a,q). We will 
discuss the form of this function for a real system and possible values of the power-law index 
f] in more detail below. 

For a neutrino-emission function of the form (!29|) . Eq. (fT3l) can be written in the form 

dN sin^dN 2qdN 2^dN 

^^^^ITB ^l^+o ri^ o -j^ = Xa"^ 6{q - qo) 30 

oR R ov oTcoiya oq Stcoi oa 



To solve this equation, we must integrate the righthand side along the characteristic. Using 
boundary condition (|2T|) . we obtain for the surface of the collapsing sphere: 



N = J Xa^+'6{q - qo)dy 



(31) 



We obtain using ([TB1)-(|251) 
II 



N 



2^coi\^Siq - qo)dq = ifcoiX-^^ 



[^[Qii - Qo] - - go]) 



Here, the function S is defined to be 



(32) 




1 



X < 0, 
X > 0. 



Note that quantities relating to the end of the characteristic are those that are measured 
by an observer located on the surface of the collapsing star. Therefore, we can omit the 
subscript "11". Substituting relations (JTSil and ( |28l) into ( |32l) . we finally obtain 



N 



3 a 



/3+I 



1 + 



(33) 



This formula specifies the neutrino distribution function measured by an observer located 
on the surface of the star and moving with this surface. 



4 PROPAGATION OF NEUTRINOS FROM THE STAR 
TO AN OBSERVER 

Let us now consider the motion of the neutrinos outside the star. The gravitational field here 
is a Schwarzschild field ^Oj; i.e., it has the metric 

The neutrino distribution function N in such a field in the spherically symmetrical case 
depends on four quantities: the angle between the neutrino trajectory and the direction 
from the center of the star 9, the neutrino energy the radius r, and the time t. We can 
find N by solving the kinetic equation, for which the boundary conditions must be obtained 
using fl551) . Let us determine the relationship between physical quantities measured on the 
surface of the star in the Friedmann and Schwarzschild reference frames. Loc al observers in 
the Schwarzschild frame are stationary, while local observers in the Friedmann frame move 
along with the stellar material. Therefore, physical quantities measured by such observers 
at some time at the same point on the stellar surface will be related by the usual Lorentz 
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transformation. The velocity in this transformation is the physical velocity of the stellar 
surface measured by a local Schwarzschild observer, which will be equal to [10] 



V 



1 
7^ 



(35) 



where h is the current Schwarzschild radius of the star. 

The circumference of the stellar surface in the Friedmann frame is 27rx, and in the 
Schwarzschild frame, 27r/i. Since dimensions perpendicular to the motion are not changed 
by the Lorentz transformation, 2ttx = 2nh, i. e., x = h. Thus, the Schwarzschild radius of 
the star is equal to x. In particular, (135|) can be rewritten as 



V 



X 



(36) 



As is shown in [S], the physical distribution function of the particles is Lorentz invariant. 
Therefore, to recalculate function (!33l) in the Schwarzschild frame, we can simply express 
the variables (g, cos??) in terms of the quantities {ws,cos9s) in the Schwarzschild fram^ 
using the Lorentz transformation and the velocity determined by (136|) . Using the results of 
[8], it is straightforward to obtain the relations: 



1 + V cos 6 s 



cost? 



Ws 



w 



y/x + COS 9s 

cos 9s + V y/x COS 9s + I 
I + V cos 9 s \fx + cos 9 s 



(37) 
(38) 



Substituting fl37j) and (l38l) into fl33l) . we finally obtain 



^0 % ' 



w 



Ix + COS 9s 



Ws - 



a/x + COS 9s 



Ws 



\Jx — 1 
a/x + COS 9^ 




-Jx COS + 1 ' 
V^ COS 9s ^ X 



X 



(39) 



Thus, we have derived the physical neutrino distribution function on the stellar surface in 
the Schwarzschild frame. T o fully specify the boundary conditions, we must determine an 
expression describing the motion of this surface. T he equation for radial free fall of a body 
in a Schwarzschild field is such that the velocity of the body is zero at infinity [10] (as was 
indicated above, the surface of the dust sphere moves in precisely this way) and has the 
form 

■v^+l\ 2 



t = In 



3 

-X2 — 



D 



(40) 



^Here and below, the subscript s denotes variables in the Schwarzschild frame specified on the surface of 
the collapsing sphere. 
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Here, D is the constant of integration. By specifying this constant, we determine the zero 
time in the Schwarzschild frame. Setting D equal to zero, we obtain 

t = \n{ ^ — - ] - -x^ -^/x (41) 



This formula describes the motion of the stellar surface. Thus, we have fully specified 
the boundary conditions for the kinetic equation in the Schwarzschild metric. The kinetic 
equation for such a field was obtained in [2], where it was shown that, in this case, the 
integrals of motion have the form 



p= sing (42) 

Vr - 1 



r 



n = w[l~-]' (43) 



^ = t- / , ^ =dz (44) 

i) (.-i)v'-^-p^(.-i) 

The first of these quantities specifies the impact parameter in the Schwarzschild field, and the 
second determines the redshift for the motion of a neutrino in such a field. The third integral 
is essentially the time when the neutrino is emitted from the stellar surface. However, it 
will be more convenient for us to use in place of the time when the neutrino is emitted, ^, 
the radius of the star at this time, x. It is obvious that these quantities are unambiguously 
related to each other in accordance with ( |4T1) : 

3 

Let us now express the quantities cos 6s and Ws in ( |39l) in terms of the integrals of motion 
Q, p and x. We substitute the values of variables corresponding to the stellar surface into 
formulas (H2l) -( l4^ (in particular, we must set r = x). This yields 



e = ln(^)-^xt-y^ (45) 



Hencfl 



P = ^^sm9s (46) 
y/x — 1 

n = ws{l-ly (47) 



Ws = J^-n (48) 
V a; — 1 



\/x^ — p'^ix — 1) 
coses = ^ (49) 



X\/X 



R . ... . . TT 

Since the surface of the coUapsing star is moving, some neutrinos emitted from the surface have ^ > 

at the time they are emitted (as a consequence of the aberration of hght). We do not take these into account 
when deriving (|49p . However, the fraction of such particles among those reaching a distant observer wiU be 
very small, so that we are quite justified in neglecting them. 
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Substituting the resulting expressions into formula 
at the stellar surface, 



for the neutrino distribution function 



N 



^0 go 



2 



(S [Q 



Here, we have introduced the notation 



u 



m 



mqo\ - 



xix 



x"^ + u 



g^mqo 



9 = 1 + 



X + u 
x"^ + u 



(50) 



(51) 



Since the space outside the star does not contain any neutrino sources, the the distribution 
function is constant along the characteristics. Therefore, the resulting expression specifies 
the neutrino distribution at an arbitrary point outside the star. 

However, we are ultimately interested in the parameters of a neutrino pulse that is 
measured by an observer infinitely distant from the collapsing star. Let the distance from 
the observer to the star he d {d ^ 1). Thus, we will be interested in the behavior of the 
integrals fH21) - fH^ for r = ci, i. e., r ^ 1. The first two integrals, ( (H2]) and fH5]) ). acquire 
the form 



p = d sin 6 

Q = w 



(52) 
(53) 



As is noted above, we use x in place of the integral ^ from system fH2|) - fl4ij) : these quantities 
are unambiguously related by fj45l) . We substitute this relation into (jHj) and express the 
time t when the neutrino reaches the radius r in terms of x and p. 



t = In 



2 3 

— T 2 

3 



2v^+ / 

i iz-l] 



1/ ■ I {z-\)^z^-p\z-\) 

We transform the integral in this formula by substituting r with d. 



dz 



(54) 



z^xlz 



i {z-\)^z^-p\z-\) 



dz = d + hi{d — 1) — a; — ln(a; — 1) + 



+ 



Zy z 



^z^-p\z-\) 



dz 



(55) 



The last integral in this expression converges as (i — > oo. Therefore, it is natural to set the 
upper limit equal to infinity. However, we can see that the initial integral of 0551) diverges. 
To elucidate the physical meaning of this fact, we substitute relation fl55l) into the initial 
formula (ISi 



t = c/ + ln(rf - 1) - 2 ln(v/x - 1) - -x^ 

3 



+ 



z- 1 



Z\ z 



z^ — p^{z — 1) 



2a/x — X + 



dz 



(56) 
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We can see that the divergence of flS^ leads to an unhmited growth in the arrival time of 
the neutrino to the observer as the distance to the observer approaches infinity.A similar 
behavior shown by function (15^ is quite natural. However, it is clear that such a dependence 
for the neutrino arrival time on the distance to the observer has absolutely no effect on the 
light curve or the observed spectrum; it only affects the arrival time of the neutrino pulse. 
Therefore, we can shift the zero time for an infinitely distant observer hj d + \n{d — 1); i.e., 
we make the substitution 



t* = t~d- \n{d - 1) 



(57) 



Substituting this value into flMl) . we obtain 

t* = -21n(Vi- 1) - 



2 5 



-X 



+ 



3 

Z\/Z 



2 — 2^/x — X + 



dz 



{51 



This expression relates the parameters t*, x and p of some phase trajectory. It can be 
considered a specification of the integral of motion x in the form of an implicit function 
x{p, t*) of the impact parameter p and the time t* of the arrival of the neutrino to a distant 
observer. Everywhere below, by x, we will mean precisely this function. Note that, since t* 
is a function only of x and p, i. e., of integrals of the motion, it is also an integral of the 
motion. 

Now, using fl52|) and fl53|) . we can readily find the number of neutrinos with energies in 
the interval [w,w + dw] passing through a sphere of radius d over a time dt*. This quantity 
can be expressed in terms of the distribution function as follows: 



$o(w,r) = STrVy pN{w,p,x{p,t*)) dp 



(59) 



where the integration is over all p values admissible for a given t* and w. Here, we have 
taken into account the equality of Q and w that follows from 0531) . Having integrated 
over the energy, we can obtain an expression for the light curve: 



Eoit*) = SttM / fmj^N{w,p,x{p,t*)) dpdw 



(60) 

We now substitute in the resulting expressions the distribution function (!50|) and transform 
them into the form 



$o(^,0 = 87r2A^ - 



PI- 

mj 



E{w ~ mqo) ~ '^{w — g'^mqo) dp (61) 




^[w — mqo ) — ^{w 



dpdw 



pm^x'^^ 2 (^2/3+11 _ 1) dp (62) 
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In all these formulas, integration over p is understood to mean integration over all values of 
p for which there exists a solution to fl58l) for a given t* . The resulting formulas specify the 
desired analytical solution for the source function of the form fl29l) . 



5 APPLICATION OF THE MODEL TO THE COL- 
LAPSE OF A REAL STAR 

We can use our model to compute the light curve and spectrum of a real collapsing star.W e 
take the mass of the star to be M = lO^M©. Since the intensity of the neutrino emission by 
the stellar matter grows with temperature, which, in turn, grows as the star is compressed, 
it is expedient to consider only the last stages of the collapse, when the vast majority of 
neutrinos are emitted. We therefore take the initial radius to be Xq = 10. 

Let us now consider the heating of the stellar matter during the compression. In our 
model, this matter is taken to be dust, which, of course, does not change its temperature 
as it is compressed. Therefore, the law describing the temperature variations of the matter 
must be specified separately. We will assume that the matter is heated like an ideal gas 
during an adiabatic compression: 

T 

= const (63) 

Here, T is the temperature of the matter and e is its density. For a relativistic gas (radiation), 
4 

7 = -. In addition, it follows from ([5]) that e oc a~^. Therefore, we can rewrite formula fl63l) 
3 

in the form 

aT = const (64) 

Denoting the initial temperature of the matter Tq (and using the fact that oq = 1), we 
finally obtain the formula 

T = — (65) 
a 

The parameters for which a star loses its stability were derived in [H], in a study of the 
evolution of massive stars. If we suppose that the subsequent compression is adiabatic and 
that the influence of the pressure can be neglected, we can use the parameters presented in 
[TT] to estimate the density and temperature when the radius of the collapsing star is equal 

to Xq. We obtain eo = 2 ■ 10^ — , Tq = 3 ■ 10^ K. 

cm 

Let us now consider the physical processes occurring up to the radiation of neutrinos 
by the heated matter. As was shown in fl2\ [13], the main mechanism for the radiation of 
neutrinos during the collapse of a supermassive star is the annihilation of electron positron 
pairs. In this case, the energy radiated per unit volume of matter per unit time is specified 
by the very simple formula ppl 

Q ~ 4 ■ IQ-'^^T^ ( (66) 
Vcm'^ s/ 



'''Expressions ([M]) . ([68]) and ([74]) are written in the cgs system. 
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Fig. 1 




-15 -10 -5 5 10 15 

Time, s 



Figure 1: Neutrino light curve. The time in the plots is the time t* for an infinitely distant 
observer whose zero time has been shifted in accordance with (1571) . For more detail about 
the choice of zero time for a distant observer, see the text following formula (l53ll . 



where the temperature is in Kelvin. Hence, /? = —9. When finding the light curve, the 
energy distribution of the neutrinos is not important and we can take all the emitted neu- 
trinos to have the same energy, go- In this case, the function /(g, a) in (|T2|) corresponding 
to the radiation intensity (|66|) and the parameters of the star presented above is given by 

2 ■ 10^^ 

f{a,q) = —a^S{q-qo) (67) 

Qo 

We thus obtain using fl62|l l^: 

Eitn=^-lO''l^il-g-')dp (^) (68) 



X 



2 



This formula specifies the desired light curve. 

Let us now consider the spectrum of the neutrinos. The spectrum of neutrinos generated 
during pair annihilations was calculated using Monte Carlo simulations in [6]; under the 



^Note that, during the calculations, all quantities in ((68|) and (j74p must be taken in the system of units 
in which c = 1 and c = 1 rg = ^^|^ = 1. The transition to cgs units is carried out via the numerical 
coefficient preceding the formula. 
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conditions of interest to us, the shape of this spectrum is approximated well by the formula 



^.2 



L{q) oc ^ (69) 

1 + exp 



i,6T ; 

where L is the energy radiated per unit volume of matter per unit time in a unit interval of 
energy. We can see that the shape of the spectrum depends on the temperature. However, 
we will assume that temperature variations affect only the total intensity of the radiation, 
leaving the spectrum unchanged; i.e., the spectrums shape will correspond to that for some 
fixed temperature Eq: 

ay 

L(q) « ^SlZ ^ (70) 

It is reasonable to choose the parameter Eq so that it is approximately equal to the tempera- 
ture of the matter when the luminosity of the star is maximum. We will take this parameter 
to be equal to the temperature of the matter when x ~ 1.43,7 g] i. e., Eq = 7To in our case. 

Thus, using expression ( l66l) for the total intensity of the radiation, we can readily show 
that the function f{q,a) in ( |T2i) is 

/(a, q) = — —a-' (71) 

' q 



q \ l + exp 



1,6^0 



Thus, in our case, the source function will be 

f{q,a)=yiq)a^ (72) 

It is straightforward to derive the form of the light curve and resulting spectrum in this 
case if we make use of the linearity of the right-hand side of f|T3|) . Expression fl72l) can be 
rewritten as _ 

/•oo 

1(1,0')= y{qo)5{q-qQ)a^dqQ (73) 







Thus, it is clear that we can obtain the spectrum corresponding to a source function of the 
form (1731) by replacing A in formula ( |6Tll with y{qo) and integrating over go- In our case, we 
obtain after straightforward manipulation: 



*(»,*•) =4. 10 r — ^^-^-^ — ^d(^^ 

.Eq 



p 



m 



1 + exp 



S(l -m^) -S(l 
w 



7. 


/ Qo 


-) 


VI, 6^0 


g m—) 
w . 


dp 



erg s 



(74) 
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Fig. 2 
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Figure 2: Spectra of the neutrinos at times t* = —10.55 s (curve with squares, luminosity 
10% of the maximum), t* = —1.25 s (curve with circles, maximum luminosity), and t* = 
7.25 s (curve with triangles, luminosity 10% of the maximum). 



The results of computations carried out using (1^ and (1711) are presented in Figs. [T] and [21 
Figure d] depicts the light curve of the collapsing star, which has a characteristic bell-shaped 

T 

appearance. The duration of the neutrino pulse is approximately 10 s, i.e., 10—. 

c 

Figure [2] presents neutrino spectra at times t* = —1.25 s (curve with circles), t* = 
— 10.55 s (curve with squares), and t* = 7.25 s (curve with triangles), with the luminosity 
being maximum in the first case and 10% of the maximum in the latter two cases. We can 
see that the spectrum becomes softer with time due to the influence of the redshift. 



6 CONCLUSION 

The model presented here can be applied with source functions ( !T2|) with fairly arbitrary 
shapes. Our approximation that the shape of the spectrum of the radiating matter does 
not depend on temperature in our derivation of formula (ITOl) was made only to simplify 
the calculation. If the source function cannot be presented in the form (1721) . it must be 
expanded in a Loran series, 

oo 
n=— oo 

and the linearity of the problem taken into account; i.e., it is necessary to obtain a solution 
for each term of the series separately and then sum the results. 
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Thus, our solution for the properties of a neutrino pulse arising during the collapse of a 
supermassive star can be used to consider various processes that generate neutrinos (plasma 
neutrinos, pair annihilation, and so forth). There is no question that the solution does not 
model a real collapse completely accurately; in particular, the assumption that the star is 
uniform during the collapse is crude. However, the resulting solution is very simple while 
also displaying all the properties that are characteristic of real systems. This makes it 
possible to obtain hght curves and spectra for such stars that are close to the real curves 
without having to carry out complex three-dimensional computations. 
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